May 97 
IFUP-TH 19/97 



Field strength correlators in full QCD [] 



M. D'Elia, A. Di Giacomo, E. Meggiolaro f\ 

Dipartimento di Fisica, 
Universita di Pisa, 
and INFN, Sezione di Pisa, 
1-56100 Pisa, Italy. 



Abstract 

We study, by numerical simulations on a lattice, the behaviour of the gauge-invariant 
two-point correlation functions of the gauge field strengths in the QCD vacuum with 
dynamical fermions. 
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1. Introduction 



A relevant role in hadron physics is played by the gauge-invariant two-point correlators 
of the field strengths in the QCD vacuum. They are defined as 

V mva (x) = (0\Tr {G, p (x)S(x,0)G ua (0)S^(x,0)}\0) , (1.1) 

where G pp = gT a G^ p is the field-strength tensor and S(x, 0) is the Schwinger phase 
operator needed to parallel-transport the tensor G ua (0) to the point x. 

They govern the effect of the gluon condensate on the level splittings in the spectrum 
of heavy QQ bound states |], |], ||. They are the basic quantities in models of stochastic 
confinement of colour j|, |j, || and in the description of high-energy hadron scattering 
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These correlators have been determined on the lattice in the quenched (i.e., pure- 
gauge) theory, with gauge group 577(2) JTl|, and also in the quenched SU(3) theory in 
the range of physical distances between 0.1 and 1 fm [[12], |13| . In this paper we compute 
them in full QCD, i.e., we also include the effects of dynamical fermions. 

The technique used is the same as in Refs. Jl2|, [HJ. The basic idea is to remove the 
effects of short-range fluctuations on large distance correlators by a local cooling procedure 



H] , T5\ . Freezing the links of QCD configurations one after the other, damps very rapidly 
the modes of short wavelength, but requires a number n of cooling steps proportional to 
the square of the distance d in lattice units to affect modes of wavelength d: 

n~kd 2 . (1.2) 

Cooling is a kind of diffusion process. If d is sufficiently large, there will be a range of 
values of n in which lattice artefacts due to short-range fluctuations have been removed, 
without touching the physics at distance d. This removal will show up as a plateau in the 
dependence of the correlators on n. 

The results are presented in Sect. 2. The determination was done at (3 = 5.35 ((3 = 
6/g 2 , where g is the coupling constant) on a 16 3 x 24 lattice with four flavours of staggered 
fermions and the Wilson action for the pure-gauge sector. We have used a standard hybrid 
Monte Carlo (HMC) algorithm, in particular the so-called $ algorithm described in detail 
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in Ref. | lq| : the trajectory length r was taken to be 0.3 with a molecular-dynamics step 
size St = 0.004. The bare quark mass was chosen to be a ■ m q = 0.01 (a being the lattice 
spacing), which should be a reasonable approximation to the chiral limit. A determination 
was also made for a ■ m q = 0.02, which we shall comment in the following. In Sect. 3 we 
discuss our results and give some concluding remarks. 



2. Computations and results 

The parametrization of the correlators is taken from Refs. [§], @, [|: 

V^ va {x) = {S^Spa - Sp a S pv ) [D(x 2 ) + £>i(x 2 ) 

, <9£>i(x 2 ) 

~\~\X pXyO p a X pXfyO p V ~\- X pX a p U X pXyO pu J — - . ^2.1) 

T> and T>\ are invariant functions of x 2 . We work in the Euclidean theory. 
It is convenient to define a T>\\(x 2 ) and a T>j_(x 2 ) as follows: 

V\\ = V + V 1 +x 2 —^ , 

ox z 

V± = V + Vx . (2.2) 

On the lattice we can define a lattice operator , which is proportional to Vp P ^ a in 

the naive continuum limit, i.e., when the lattice spacing a — > 0, |13|1 . Making use of 
the definition (2.2) we can thus write, in the same limit, 

Vf(da) ~ a 4 P||(rfV) + C(a 6 ) , 

V L Ada) ~ a 4 P ± (rf 2 a 2 ) + C(a 6 ) . (2.3) 

Higher orders in a in Eq. (2.3) as well as possible multiplicative renormalizations are re- 
moved by cooling the quantum fluctuations at the scale of the lattice spacing, as explained 
in the Introduction. 

The only scale in our system is the lattice spacing a: its value in physical units depends 
on f3. Here we work with only one value of /3, so we could present our results directly in 
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units of a. However, in order to facilitate the comparison with our previous works [12, 13 
we shall use the familiar parametrization 



a(P) = 4-fiP) , (2-4) 
where the scaling function f(/3) is given by the usual two-loop expression: 

, o x 231/625 , a v 

for gauge group SU(3) and Nf = 4 flavours of quarks. in Eq. (2.4) is an effective 
A-parameter for QCD in the lattice renormalization scheme, with Nf = 4 flavours of 
quarks. With this parametrization: 



Pi/I/?)"* = il?/ ( ' 3) ) ■ (2 ' 6) 

We have measured the correlations on a 16 3 x 24 lattice at distances d ranging from 3 to 
8 lattice spacings and at (3 — 5.35. At this value of (3 the lattice spacing a(/3), extracted 



from the string tension or the p mass, is of the order of 0.11 fm ||T7| , [TSf] , so that the 
lattice size is approximately 2 fm and therefore safe from infrared artefacts. In Fig. 1 we 
display the results for Vff(f3)~ 4 and V^f((3)~ 4 versus G? p h ys = (d/Ap) f((3), for a quark 
mass a ■ m q = 0.01. Measurements have been done on a sample of 150 configurations, 
each separated by 15 HMC trajectories. Statistical errors have been estimated by using 
a standard blocking procedure. As in Ref. |T3| we have tried a best fit to these data with 
the functions 

V(x 2 ) = A exp(-\x\/X A ) + -^-exp (-|s|/A a ) , 



kr 



X?!(x 2 ) = A 1 exp(-|x|/A A ) + ^Texp(-|x|/A a ) . (2.7) 



We have obtained the following results: 



^ = (1.74 ± 0.24) x 10 10 , ^ = (0.20 ± 0.10) x 10 10 



a = 0.71 ±0.03 , ai = 0.45 ±0.03 , 

544 ±27 , —— = 42 ±11, 2.8 



Xa^f A a Aj 
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with x 2 /N<i.o.f. — 0.5. The continuum lines in Fig. 1 correspond to the central values of 
this best fit. 

The corresponding results for the quark mass a ■ m q = 0.02 are displayed in Fig. 2, 
using for Ap the same value adopted for a ■ m q = 0.01: we shall comment on this point in 
the next Section. Measurements have been done on a sample of 30 configurations, each 
separated by 20 HMC trajectories. A best fit to these data with the same functions (2.7) 
gives the results 

^ = (3.48 ± 0.42) x 10 10 , ^ = (0.46 ± 0.21) x 10 10 , 

a = 0.66 ±0.03 , ai = 0.39 ± 0.03 , 

— 1— = 631 ±23 , — i- = 61 ±20, (2.9) 
XaAf X a A F 

with x 2 / N d.o.f. ^0.7. A gain, the continuum lines in Fig. 2 correspond to the central 
values of this best fit. 



3. Discussion 

Two quantities of physical interest can be extracted from our lattice determinations: 

1) the correlation length A^ of the gluon field strengths, defined in Eq. (2.7); 

2) the so-called gluon condensate, defined as 

G * = : G%G% :) = (3- 1 ) 

Both of them play an important role in phenomenology. The correlation length is 
relevant for the description of vacuum models 0, [5], |(|. The relevance of the gluon con- 
densate was first pointed out by Shifman, Vainshtein and Zakharov (SVZ) |]l9fl . It is a 
fundamental quantity in QCD, in the context of the SVZ sum rules. 

From lattice we extract Xa in units of lattice spacing a. To convert these units to 
physical units, the scale must be set by comparison with some physical quantity. This 
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was done in Refs. [jlT], [H| by computing the string tension and the p mass on the lattice 
and comparing them with the physical values. Their lattice was identical to ours (16 3 x 24) : 
they also used the same value of j3 (5.35), as well as Nf = 4 flavours of staggered fermions 
with mass a-m q = 0.01, as we did. Their estimate for the lattice spacing is a ~ 0.11±0.01 
fm. This gives: 

X A = 0.34 ± 0.02 ± 0.03 fm (a • m q = 0.01) . (3.2) 

The first error comes from our determination [Eq. (2.8)], the second from the error in 
converting the lattice spacing to physical units. From a ■ m q = 0.01 to a ■ m q = 0.02 
the value of the effective Ap, defined by Eqs. (2.4) and (2.5), can change in principle. 
Anyway, no published determination of A^ (i.e., of the lattice spacing in physical units) 
exists for a ■ m q = 0.02. Some data on the pseudoscalar and vector meson masses, for 



quark masses a ■ m q larger than 0.01, have been published in Ref. [|20|. We have tried 



to extract A^ (i.e., the lattice spacing) from those data by using the same procedure of 



Ref. \ T% . We estimate that from a ■ m q = 0.01 to a ■ m q = 0.02 the effective mass-scale 
Ap does not change appreciably within the errors. Assuming, as an indication, the same 
value of Ap as for a ■ m q = 0.01, we then get: 

\ A = 0.29 ± 0.01 ± 0.03 fm (a ■ m q = 0.02) . (3.3) 

The values (3.2) and (3.3) must be compared with the quenched value |1^, 

X A = 0.22 ± 0.01 ± 0.02 fm (YM theory) . (3.4) 

Here the value A^ ~ 4.9 ± 0.5 MeV has been assumed for the pure-gauge A-parameter 



21| . The correlation length increases when going from chiral to quenched QCD and 
this tendency is confirmed by the fact that decreases by increasing the quark mass. 
Of course, a precise determination of A^ should be done with more realistic values for the 
quark masses. 

We now come to the gluon condensate. Our lattice provides us with a regularized de- 
termination of the correlators. At small distances x a Wilson operator-product- expansion 
(OPE) |22| is expected to hold. The regularized correlators will then mix to the iden- 
tity operator 1, to the renormalized local operators of dimension four, 9L ^'-G a VLV G° VLV : and 
mf.qfqf. (/ = 1, . . . , Nf, Nf being the number of quark flavours), and to operators of 
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higher dimension: 



N r 



;V„(x) ~ C 1 (x){l) + C g {x)G 2 + £ C f (x)m f (: q f q f :) + ... . (3.5) 



2tt 2 ^ y ' x^o { , 

The mixing to the identity operator Ci(x) shows up as a c/|x| 4 behaviour at small x. The 
mixings to the operators of dimension four C g (x) and Cf(x) are expected to behave as 
constants for x — > 0, while the other Wilson coefficients in the OPE (3.5) are expected to 
vanish when x — > (for dimensional reasons). The coefficients of the Wilson expansion 
are usually determined in perturbation theory and are known to be plagued by the so- 
called infrared renormalons (see for example Ref . []23| and references therein) . In our case 



this means that, due to the infrared renormalon pole, terms coming from the mixing to 
the identity operator can produce by resummation a term which simulates a mixing to a 
condensate | 23| . There is no specific recipe to disentangle this u perturbative 11 contribution 



to the condensates from a possible "genuine" value of them. A similar problem is present 
in any Wilson OPE in QCD, and in particular in the expansion which leads to the SVZ 
sum rules [fL9] . 

A practical way out, which provides good results for the sum rules, is to assume that 
the leading perturbative determination of the Wilson coefficients is a good approxima- 
tion to the unknown determination which should be done by perturbing around the real 
vacuum of the theory (see for example Ref. [23| for a discussion about this point). In 



this spirit, we shall assume that the renormalon ambiguity can be safely neglected in the 
extrapolation for x — > of our correlators. With the normalization of Eq. (3.5), this 
gives C g (0) ~ 1. On the same line, the contribution from the quark operators in (3.5) can 
be neglected, because the corresponding condensates q/qf :) are much smaller than 
G 2 and the mixing coefficients Cf(x) are of higher order than C g (x) in the perturbative 
expansion. Within these approximations, we get the following expression for the gluon 
condensate, in terms of the parameters defined in Eq. (2.7): 

G 2 --^-(A + A 1 ) . (3.6) 
At a ■ m q = 0.01 this gives, in physical units, 

G 2 = 0.015 ± 0.003 tomt Gey4 (a-m q = 0.01) . (3.7) 
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At a ■ m q = 0.02 we obtain: 



G 2 = 0.031 ± 0.005 t°0M7 GeV 4 (a ■ m g = 0.02) . (3.8) 

The two values (3.7) and (3.8) should be compared with the corresponding value in the 
quenched theory [Tj| : 

G 2 = 0.14 ± 0.02 tom Gey4 ( YM theory) . (3.9) 

The gluon condensate G 2 appears to increase with the quark mass, as expected, tending 
towards the asymptotic (pure-gauge) value of Eq. (3.9). Contrary to the previous dis- 
cussion for Xa, we have here a theoretical tool to understand the dependence of G 2 on 
the quark masses. According to Ref . |2f| , we expect the following low-energy theorem to 
hold for small quark masses (m/ -C fi, being the renormalization scale): 

^ : G % G % :) = -|<: ~ Wt :) , (3.10) 
where b = 11 — l-ATy, for a gauge group SU(3) and Nf quark flavours. For the two 



gluon condensates (3.7) and (3.8) one must use the renormalized quark masses uif [26 
corresponding to a ■ m q = 0.01 and a ■ m q = 0.02 respectively: for a ■ m q = 0.01 we have 
approximately m f ~ 44 Me V. Making use of the popular values for the quark condensate 
((qq) — —0.013 GeV 4 ^|) and for the physical quark masses (m u ~ 4 MeV, m d ~ 7 
MeV and m s ~ 150 MeV), we can thus extrapolate from the value (3.7) to the physical 
gluon condensate, obtaining: 

G^ hysical) ~ 0.022 GeV 4 . (3.11) 



The procedure used is the same as in Ref. [25|. The prediction (3.11) agrees with the 
empiric value obtained from experiments E7\, E8[: Q^ m P mc ^ ~ 0.024 ± 0.011 GeV 4 . 
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FIGURE CAPTIONS 



Fig. 1. The functions T>±f(f3) 4 (upper curve) and T>^f(f3) 4 (lower curve) versus physical 
distance, for quark mass a ■ m q = 0.01. The curves correspond to our best fit [Eqs. 
(2.7) and (2.8)]. 

Fig. 2. The same as in Fig. 1 for quark mass a ■ m q = 0.02. The curves correspond to our 
best fit [Eqs. (2.7) and (2.9)]. 
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